Multi-band spectroscopy of ultracold fermions: Observation of reduced tunneling in 

attractive Bose-Fermi mixtures 
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We perform a detailed experimental study of the band excitations and tunneling properties of 
ultracold fermions in optical lattices. Employing a novel multi-band spectroscopy for fermionic 
atoms we can measure the full band structure and tunneling energy with high accuracy. In an 
attractive Bose-Fermi mixture we observe a significant reduction of the fermionic tunneling energy, 
which depends on the relative atom numbers. We attribute this to an interaction-induced increase 
of the lattice depth due to self-trapping of the atoms. 

PACS numbers: 37.10.Jk, 67.85.Pq, 71.10.Fd, 78.47.-p 



Quantum gases in optical lattices have proven to be 
a versatile model system to investigate unsolved many- 
body problems from condensed-matter theories in a 
highly controlled fashion [T]. The tunability of interac- 
tion strength, lattice depth, and lattice geometry makes 
these systems particularly well suited to explore depen- 
dencies and limits of anticipated models like the Hub- 
bard Hamiltonian P] [3] . Especially fermionic quantum 
gases have attracted much interest over the last years due 
to their close analogy to electrons in crystals. This led, 
e.g., to the observation of Fermi-surfaces j^, antibunch- 
ing [5], and the fermionic Mott insulator [6l[7]. At the 
same time quantum gases in optical lattices allow to real- 
ize completely new systems like mixtures of bosonic and 
fermionic atoms where several exotic phases have been 
predicted. Among them are density waves f5^, supersolids 
[9], and polaronic quasiparticles ^lOj . Pioneering experi- 
ments found a pronounced shift of the bosonic Mott in- 
sulator transition |11H13| and more recently a renormal- 
ization of the on-site interaction [T3]. The realization 
of these systems over a broad range of parameters has 
extended the theoretical interest. There is a controver- 
sial debate on the correct theoretical description of Bose- 
Fermi mixtures: Various models have been discussed, in- 
cluding adiabatic heating during the lattice loading [151 
and an effective potential approach including the admix- 
ture of higher bands |13| ITS] . A recent ansatz introduces 
an extended Hubbard model with nearest neighbor inter- 
action 117] . For a complete comparison with experiments, 
it is important to independently measure both the effec- 
tive on-site interaction and the effective tunneling in the 
mixture, where the latter has not been directly observed 
in experiments yet. 

In this Letter, we study the tunneling properties of ul- 
tracold fermions in optical lattices with high accuracy. 
This is done for pure spin-polarized fermions and for at- 
tractive mixtures of fermionic Potassium and bosonic Ru- 
bidium atoms. As a key result we directly observe an 
interaction-induced change in the effective lattice depth 
and thus a reduction of the fermionic tunneling energy J 
due to the bosonic component [TE] . We explain this in an 



effective potential picture where mutual self-trapping of 
both species is expected |13lfTB] . For these measurements 
we apply a novel multi-band spectroscopy for fermions in 
optical lattices based on lattice modulation and a subse- 
quent band mapping. Similar to ARPES in solid state 
physics, our method can accurately resolve the excitation 
spectrum with full momentum resolution. This allows for 
the determination of the lattice depth and the tunneling 
energy with high accuracy but is also promising to study 
excitons as well as non-equilibrium dynamics of particle 
and hole excitations. 

We prepare a spin-polarized ultracold mixture of 
bosonic ®^Rb and fermionic ^°K atoms in a crossed op- 
tical dipole trap with typically 10^ fermions at a tem- 




Figure 1. (color online), (a) Experimental sequence, (b) 
Experimental procedure in momentum space, (c) Typical 
absorption image and column sum with nonresonant mod- 
ulation. The atoms occupy the first Brillouin zone, (d) With 
resonant modulation: Particle-hole excitations are clearly vis- 
ible. 
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Figure 2. (color online). Momentum-resolved band struc- 
ture of non-interacting fermions in an optical lattice of 5 E,- : 
Shown are the column densities of the momentum distribu- 
tions for different modulation frequencies. Central plateau 
are atoms in the first Brillouin zone. Missing particles are 
holes, representing the reduced zone scheme. Narrow peaks 
at higher momenta are the outcoupled particles, representing 
the extended zone scheme. 



perature below 0.2 Tp and 2 x 10^ atoms in the BEC. 
The dipole trap is operated near the magic wavelength at 
811 nm, compensating the differential gravitational sag, 
with an average trap frequency of lii = 27r x 50 Hz. We 
superimpose a 3d optical lattice at A = 1030 nm which we 
ramp up in 100 ms. Note that the Rubidium atoms expe- 
rience an approximately 2.5 times stronger lattice than 
the Potassium atoms with respect to their recoil energies 
due to the different masses and detunings. In the fol- 
lowing the lattice strength will be denoted in fermionic 
recoil energies Ej- — h^k^/2ra with the mass m of '^''K 
and the lattice vector k]^ = 27r/A. Depending on the 
lattice depth, the fermionic ground state is a metal or 
band insulator [H [7] while the bosons form a superfluid 
or Mott insulator. 

The spectroscopy is performed as follows (see Fig- 
ure IT]) : We excite the system to higher bands by modu- 
lating the depth of the optical lattice in one direction for 
1 ms with a variable frequency and a typical modulation 
depth of 15% [ini Uni- Due to the different curvature 
of the individual bands, only a specific quasimomentum 
class of atoms can be resonantly excited at a fixed mod- 
ulation frequency. Finally we map the quasimomentum 
distribution onto real momenta by ramping down the op- 
tical lattice nonadiabatically in 200 jis suppressing band 
transitions [21] . The ramp is designed to be much faster 



than the trap dynamics to prevent redistribution within 
each band. Due to the finite ramp time, however, the 
distribution is smoothed out at the zone boundaries [22] 
[see Figure llTc)]. We investigate the resulting momen- 
tum distribution after a time-of-flight of typically 20 ms 
thus spatially separating different bands and momenta. 
Applying off-resonant modulation, the atoms stay in the 
ground state and occupy the first Brillouin zone [see Fig- 
ure [l[c)[. If the modulation frequency corresponds to a 
transition energy, a specific momentum class is excited 
to a higher band and mapped onto its corresponding 
Brillouin zone creating a hole within the first Brillouin 
zone [see Figure [lTd)[. Note that in contrast to local- 
ized particle-hole states in the Mott insulating regime |6] 
we create particle-hole excitations in momentum space. 
Varying the modulation frequency gives access to the full 
band structure with full momentum resolution. 

To demonstrate the capabilities of this method, we first 
apply it to the well-known case of a pure non-interacting 
fermionic sample in an optical lattice as shown in Fig- 
ure [2] Experimentally, we measure the energy difference 
between the first band and the excited bands, obtaining 
all information to deduce the full band structure. The 
outcoupled particles correspond to the extended zone 
scheme of the band structure. The rise of the band 
gap and the corresponding flattening of the dispersion 
are clearly visible at the edge of the Brillouin zones. 
The holes represent the quasimomenta of the outcoupled 
atoms and thus the band structure within the reduced 
zone scheme in a textbook-like manner. Note the in- 
verted curvature of the second band, which becomes ap- 
parent in the reversed dispersion of the holes. To our 
knowledge, this is the first experimental observation of 
the full band structure for ultracold fermions in optical 
lattices. 

To analyze the band structure quantitatively, we deter- 
mine the center of mass of the outcoupled atoms in the 
time-of-flight pictures and extract the corresponding mo- 
menta for all applied modulation frequencies. The result 
is shown exemplarily in Figure Wia) for the third band 
at a lattice depth of llEi- with the fermions forming a 
pure band insulator. The momentum transfer is given 
in units of kez = 27r/A. We estimate the lattice depth 
of the system by a least square fit of a homogeneous Id 
band structure (red solid line) to the experimental data. 
In our data, we clearly observe a deviation of the exci- 
tation spectrum from the homogeneous case particularly 
at small and large momenta. We find good qualitative 
agreement of this observation comparing our data to a 
numerical linear response calculation of a finite Id sys- 
tem including the harmonic confinement. Especially the 
bending of the spectroscopic signal at the zone edges is 
well reproduced. In an intuitive picture we interpret this 
effect as the coupling of different momentum components 
due to the inhomogeneity especially in the first band. 
This coupling influences the eigenstates and the density 
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Figure 3. (color online), (a) Extracted dispersion of the third 
band. Red circles are the centers of mass of the excited atoms. 
The dashed line is an inhomogeneous Id calculation, repro- 
ducing the data qualitatively. The red solid line is a fit to 
the data for 0.55 < q < 0.85 kez with a homogeneous Id 
band structure. Shaded areas show the dispersion for slightly 
different lattice depths, (b) Histogram showing the relative 
number of outcoupled atoms. 



of states at the upper and lower edge of the first band. 
In these parts of the spectrum no atoms are present to 
be transferred into the third band. For a specific range 
of momenta the effect of the inhomogeneity is negligible. 
Therefore we restrict the fitting to our data to a momen- 
tum range between 0.55 and 0.85 kez- From this fit we 
can determine the lattice depth with very high accuracy. 
The dominant errors of about 2% are of systematic na- 
ture, given by the spatial dependence of the lattice depth 
due to the gaussian shape of our lattice beams. For non- 
interacting particles, the lattice depth fully determines 
the tunneling energy. The fitted value of 11.3(2) Ej- cor- 
responds to J/h — 67(3) Hz. In addition to the dispersion 
relation we count the number of transferred atoms in the 
corresponding Brillouin zone, shown as a histogram in 
Figure p[b) . It has a maximum at low momenta and is 
small where the signal bends away from the homogeneous 
dispersion as expected. 

We now turn to our main results, investigating the in- 
fluence of interspecies interaction on the fermionic band 
structure: We apply our spectroscopy method to study 
an attractive mixture of bosonic *^Rb atoms in the state 
|1, 1) and fermionic ^°K atoms in the state |9/2,9/2). In 
this system, interaction effects represented by the strong 
attractive interspecies scattering length abf ~ — 215ao 
between *°K and ^^Rb [53] and the repulsive bosonic 
scattering length abb ~ 100 ao, play a crucial role. We 
choose lattice depths such that the Rubidium atoms form 
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Figure 4. (color online), (a) Sketch of the effective potential 
for different bosonic fillings. Fermions are more tightly bound 
if attractive bosons are added, (b) Shift of the outcoupled 
momentum against the relative atom number N-h^/Nk at 5 Er 
and fixed modulation frequency of 34 kHz. Lower outcoupled 
momenta correspond to a deeper lattice. The shaded area is 
a guide to the eye. (c) Fermionic band structure for different 
bosonic admixtures. Solid lines are fitted dispersions, (d) 
Histogram showing the relative number of outcoupled atoms. 



a Mott insulator and the Potassium atoms a band insula- 
tor. Due to the 2.5 times higher lattice depth the tunnel- 
ing energy of the bosons is very small, while the fermions 
experience a moderate lattice potential with high tun- 
neling energies. We describe this situation in an effective 
potential picture |13l I16| : If a boson and a fermion are 
at the same lattice site, the attractive interaction leads 
to an enhanced localization and thus a suppression of the 
tunneling. This can be described as an effective poten- 
tial which is proportional to the respective density of the 
other component. Since both species mutually influence 
each other, this is called self-trapping. The local defor- 
mation of the wavefunctions can be interpreted as the 
admixture of higher bands. The resulting potential for 
the fermionic atoms due to the admixture of bosons is 
sketched in Figure H[ a). In a homogeneous system with 



a fixed number of Rubidium atoms per site the effective 
potential for the fermions has the same periodicity as 
the original lattice but the effective lattice depth is sig- 
nificantly increased. Here, for the first time we are able 
to directly observe this interaction shift for different mix- 
ture ratios of Potassium and Rubidium. To maximize the 
signal-to-noise ratio, we restrict our measurements to the 
third band, which has the highest excitation rate. 

A detailed scan of the third band is shown in Fig- 
ure BJc). The purely fermionic system (red circles) at 
a lattice depth of 7.6 E^ shows the typical characteristics 
at low momenta, that can be qualitatively understood as 
in Figure [3] The admixture of bosons leads to a clear 
shift of the band structure and a suppression of the ex- 
citation fraction. For moderate atom numbers with a 
bosonic Mott insulator of mainly unity filling, there is a 
small but distinct overall shift of the spectroscopic sig- 
nal (green diamonds) . Using the same fit as for the pure 
fermionic case, we obtain a lattice depth of 8, 2Ei.. This 
is consistent with the theoretical prediction without cor- 
relation effects. For higher particle numbers a shell struc- 
ture for the bosons emerges. Here we estimate a maximal 
occupation of two in the Mott insulator. In this case, 
the spectroscopic signal shows large deviations from the 
simple model (blue squares). We attribute this to the 
increased inhomogeneity of the effective potential due to 
the bosonic shell structure. If we fit a band structure as 
before we obtain a lattice depth of 9.2 Ei- which is signif- 
icantly higher than the uncorrelated shift for two bosons 
per site. Assuming the fermions as still non-interacting, 
we can derive their tunneling energy from the lattice 
depth with a simple calculation. We obtain values for 
J/h of 160 Hz, 138 Hz, and 109 Hz for the three respec- 
tive cases showing a clear reduction of the fermionic tun- 
neling due to the admixture of bosons. The number of 
excited atoms is also consistent with our interpretation: 
Since the fermions experience a higher lattice depth in 
the presence of bosons, less atoms are transferred at low 
energies, as clearly visible in Figure llFd). 

To investigate the particle-number dependence in more 
detail, we performed a measurement at fixed modulation 
frequency for a shallow lattice of 5Er [see Figure l4[b)]. 
For pure fermions, the modulation frequency of 34 kHz 
results in a momentum of fc = 0.58 kez- For an increas- 
ing admixture of bosons the mean outcoupled momentum 
decreases. This is consistent with the effective potential 
approach which predicts an increase of the effective lat- 
tice depth as in Figure l4[c) . The non- linear slope indi- 
cates enhanced localization due to stronger correlation 
effects for higher particle numbers. 

In conclusion, we have performed a fully momentum 
resolved spectroscopy of ultracold fermions in optical lat- 
tices for the first time. We have investigated a mixture 
of fermions and localized bosons in an optical lattice and 
observed an interaction induced shift of the fermionic 
band structure due to the attractive bosons as expected 



from an effective potential approach. The effective lattice 
depth was increased up to 20% while the corresponding 
tunneling energy was reduced by more than 30%. Further 
we found evidence for a correlation induced enhancement 
of the interaction shift. As a complementary measur- 
ment to |14| . this is an important step for the compari- 
son with theory. In agreement with |13L I14| , we conclude 
that a simple single-band Hubbard model is not sufficient 
for the realistic description of strongly attractive quan- 
tum gas mixtures. Moreover, we observe a confinement- 
induced deformation of the spectroscopy signal at small 
momenta, consistent with numerical calculations. Our 
easy-to-implement spectroscopy method is also promising 
for the investigation of interacting fermi gases in optical 
lattices. It is in particular suited for band gap excitations 
such as excitons and the dynamics of particle and hole 
excitations. 
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